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The phase transition from hadronic to quark matter may take place already during the early 
post-bounce stage of core collapse supernovae when matter is still hot and lepton rich. If the phase 
transition is of first order and exhibits a barrier, the formation of the new phase occurs via the 
nucleation of droplets. We investigate the thermal nucleation of a quark phase in supernova matter 
and calculate its rate for a wide range of physical parameters. We show that the formation of the 
first droplet of a quark phase might be very fast and therefore the phase transition to quark matter 
could play an important role in the mechanism and dynamics of supernova explosions. 
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I. INTRODUCTION 

The possibility that a first-order phase transition in 
dense matter has implications for the explosions of su- 
pernovae was first proposed by Migdal et al. about 30 
years ago Since then, a large number of papers ad- 
dressed this issue with different hypotheses on the na- 
ture of the phase transition: pion or kaon condensate 
matter or quark matter, different models to compute the 
equation of state (EoS), different simplifications for the 
complex hydrodynamical evolution of collapsing massive 
stars and different approaches for the neutrino Boltz- 
mann transport However, only very recently was it 

possible to perform simulations using general relativistic 
Boltzmann neutrino transport equations and adopting 
realistic equations of state for quark matter 0]. In 
Ref. 0, in particular, it was shown that the phase tran- 
sition to quark matter can occur already during the early 
post-bounce phase of a core collapse supernova event and 
that it produces a second shock wave (the first being the 
usual shock wave after the bounce) which triggers a de- 
layed supernova explosion, even within spherical symme- 
try, for masses of the progenitor star up to 15 Mq. The 
formation of quark matter is also responsible for the emis- 
sion of a neutrino burst, typically a few hundred millisec- 
onds after the first neutronization burst, which could be 
detected by currently available neutrino detectors, repre- 
senting a spectacular possible signature of quark matter 
formation in compact stars (see also 01 )• 

In the studies mentioned above, an important physical 
phenomenon is neglected for the sake of simplicity: the 
process of phase conversion in a first-order transition is 
actually driven by the nucleation of finite-size structures, 
such as droplets or bubbles, of the new phase within the 
old phase. The surface tension, cr, is the physical quantity 
that determines the nature of the process of phase con- 
version. If a is sufficiently small, nucleation can be very 
fast and the new phase is produced almost in mechani- 
cal, thermal and chemical equilibrium with the nuclear 



matter phase. An intermediate value of a might render 
nucleation very difficult and the nuclear phase can be 
metastable for a significant amount of time. Then, the 
process of formation of the new phase, once triggered, 
would be a genuine nonequilibrium process, in which dif- 
ferent mechanisms can take place: deflagration, detona- 
tion, convective instabilities (see, e.g., Refs. and ref- 
erences therein) . Finally, nucleation is highly suppressed 
for large values of cr, and the formation of the new phase 
may only proceed via spinodal decomposition, if the den- 
sity achieved is high enough to flatten out the activation 
barrier. (See Ref. [3 

for a recent detailed discussion of 
the phase conversion process in a first-order phase tran- 
sition.) 

The nucleation of quark matter in neutron stars has 
been explored mainly within a scenario, proposed in 
Ref. ll| , in which the formation of quark matter occurs 
only when the protoneutron star is almost completely 
deleptonized and the temperature has already dropped 
to, say, 1 MeV. Under these conditions, quantum nucle- 
ation has been shown to be the most important mech- 
anism for the formation of a quark phase [l3 - il6j . also 
when color superconducting quark phases are present 
pTi . [Tsj . A less explored scenario is the nucleation of 
quark matter in hot and lepton-rich protoneutron stars. 
Refs. [l^, HO] present the first estimates and calculations 
showing thermal nucleation to be very efficient for tem- 
peratures of roughly 10 MeV and practically negligible 
for temperatures below 2 MeV. The simplifying assump- 
tion adopted in those papers (no Icptons are included in 
the quark equation of state) might be, however, not real- 
istic considering that for the large temperatures needed 
to nucleate quark matter the neutrino mean free path 
is small, and therefore neutrinos are trapped. In the 
presence of neutrinos, the critical densities for the phase 
transition are shifted toward larger values compared to a 
deleptonized neutron star. Moreover, as we will discuss 
in the following, the assumption of flavor conservation 
during the phase transition (also adopted in Ref. (2l| ) 
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might be too conservative in light of the quark density 
fluctuations that are evidently present. 

The main goal of this paper is to compare the time 
scale associated with the phase conversion with the dy- 
namical time scale of core collapse supcrnovae during 
which quark matter might be eventually formed. In par- 
ticular, since the explosion process occurs within a time 
scale of few hundred ms, the nucleation time must be of 
the same order of magnitude if quark matter plays indeed 
a role in the explosion mechanism. 

For this purpose, we reconsider thermal nucleation of 
quark matter within the scenario proposed in Ref. of a 
phase transition occurring already during the early post- 
bounce stage of a core collapse supernova. This implies a 
value for the critical density < 2nQ (where rig = 0.16 
fm"'^ is the nuclear matter saturation density). We sys- 
tematically investigate the windows of free parameters of 
the model adopted to compute the equations of state and 
the corresponding nucleation time scale. Our strategy 
is to present underestimates of this time scale, so that 
a successful phase conversion could somehow constrain 
the equation of state parameter space. We also discuss 
the important issue of flavor number conservation during 
the phase transition. We argue that thermal nucleation 
might indeed be efficient under the conditions realized in 
a star soon after bounce. In such case, the phase transi- 
tion would proceed very fast, thus confirming the results 
found in Ref. Q. 

The paper is organized as follows. In Sec. II we present 
the model we adopt for the equation of state, as well as 
a brief self-contained description of homogeneous nucle- 
ation and the role of statistical fluctuations. In Sec. Ill 
we present our results for the nucleation of nonstrange 
and strange quark matter. There, we also discuss quan- 
tum nucleation and the spinodal instability. Besides, we 
verify our equations of state by computing the stellar 
structure that emerge form the Tolman-Oppcnheimcr- 
Volkov (TOV) equations. Finally, Sec. IV contains our 
conclusions. 



II. PHENOMENOLOGICAL FRAMEWORK 
A. Equations of state 

The phase transition from nuclear to quark matter is 
implemented, as customary, by matching the equations 
of state for each phase. In order to do that, one has to 
choose appropriate models to compute the equation of 
state for each of the two phases and impose the condi- 
tions for mechanical, thermal and chemical equilibrium 
to determine the transition point. For nuclear matter, 
we adopt the relativistic mean field model with the TMl 
parametrization p^ . often used in supcrnovae simula- 
tions. For quark matter, we choose the MIT ba g rn odel 
[23I [23 | including perturbative QCD corrections |25l [26j . 

Moreover, we consider two types of EoS for quark mat- 
ter, both including the pressure from electrons and neu- 



trinos, which are still present at this point of the stel- 
lar evolution. The first EoS contains only up and down 
quarks, while in the second we include a massive s quark 
as well. The quark model has as free parameters the 
bag constant -B, the mass of the strange quark mg (when 
present) , and the value of the coefficient c, that accounts 
for the perturbative QCD corrections to the free gas pres- 
sure as follows: 
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where B is the bag constant, p,, is the contribution from 
the (massive) strange quark 
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and terms 0{m'^/^'^) ^ 1% in Eq. ([2]) have been ne- 
glected. Notice that, for the u-d equation of state, we 
neglect the terms related to the strange quark. 

The free parameters are fixed by requiring a critical 
density for the phase transition below two times the nu- 
clear saturation density ng = 0.16 fm~'^ for the typical 
conditions of matter in the core of a star during a su- 
pernova collapse, i.e. temperatures T = 10 — 20 McV 
and lepton fractions = 0.3 — 0.4. In this way, we ful- 
fill our initial hypothesis of formation of quark matter in 
the early post-bounce stage. Another important crite- 
rion to fix our free parameters comes from the computa- 
tion of the maximum mass of cold hybrid stars. Taking 
into account the recent measurement of the mass of PSR 
J1903-h0327, M = (1.671 ± O.OO8)M0 we require a 
maximum mass in agreement with this value. Finally, 
we investigate two possible scenarios for the appearance 
of strange quarks in the system. Since the nuclear EoS 
does not contain strangeness, a phase transition directly 
to strange quark matter might be difficult if wc consider 
the slowness of weak reactions producing strange quarks 
with respect to the fast deconfinement /chiral phase tran- 
sition process driven by the strong interaction. There- 
fore, we discuss a first case in which the phase transition 
involves two-flavor quark matter (strange quarks will be 
produced only later, via weak interaction, as suggested 
in Ref. [1^). In the second scenario, we consider a fast 
production of strange quarks: since we assume critical 
densities of the order of 2 times the saturation density 
and temperatures of a few tens of MeV, it is possible 
that a small seed of strange matter appears in the sys- 
tem through hypcrons or kaons [2^ . Once strangeness is 
produced in the hadronic matter, this would trigger the 
phase transition directly to strange quark matter. 

The equations of state for nuclear matter and quark 
matter are calculated under conditions of local charge 
neutrality, local lepton fraction conservation (i.e., the two 
phases have the same 1^), and weak equilibrium. Under 
these assumptions, the conditions of phase coexistence, 
as found in Ref. [30| , are the equality of the total pressure 
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of the two phases — P'^ and the following condition 
of chemical equilibrium: 



flu + 2 fid + = ^^eff, 



(3) 



where fin and fj,^ are the chemical potentials of neu- 
tron and neutrinos within the nuclear phase, and fi^, 
fid and fi^ are the chemical potentials of up and down 
quarks and of neutrinos within the quark phase, respec- 
tively. The quantity fif,f f is an effective chemical poten- 
tial, which is always the same in both phases. Notice that 
the condition ([3]) is always valid, although the condition 
P^ = P^ is valid only in the transition point. Here we 
use the zero-temperature equations of state, since a tem- 
perature of the order of a few tens of MeV does not alter 
considerably the equation of state ^ . 

Moreover, we assume that the different degrees of free- 
dom in both phases are in chemical equilibrium with re- 
spect to weak reactions (see Sec. II C): 
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Finally, the conditions of local charge neutrality and 
local lepton fraction within the two phases allow us to 
compute all chemical potentials in terms of only one in- 
dependent chemical potential: 
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where n,; (i = n,p,u,d^s^e^v) are the densities of the 
different species of particles. 

Notice that fixing locally results in a jump of the 
chemical potential of neutrinos at the interface of the 
phase transition. This might be not very realistic, as dis- 
cussed in Ref. 
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where it has been shown that a mixed 
phase should instead be considered due to the global con- 
servation of the lepton number. We leave the discussion 
of nucleation in mixed phase for a future study. 



B. Thermal homogeneous nucleation 

First-order phase transitions are very well known even 
from simple everyday examples, such as in the melting of 
ice. In such a case, the conversion from one phase to the 
other usually occurs slowly and very close to the thermo- 
dynamical equilibrium, following the so-called Maxwell 



^ For a free massless gas, the corrections would be 0{T'^ / ^j?) 
1%. 



construction. However, when some relevant external con- 
trol parameter (such as the temperature or the density) 
changes abruptly when a system is near the transition, 
the system finds itself in an unstable situation. For def- 
initeness, consider a system initially homogeneous in a 
low-density phase (a "gas"), and close to the transition 
line to a high-density ("liquid") phase. Now, let it suf- 
fer a sudden compression. Although the system was pre- 
pared at the gas phase, the free energy at the new, higher 
density disfavors the gas phase and the liquid phase now 
becomes the stable one: phase conversion is about to 
begin. 

The ever-present thermal and quantum fluctuations 
will not be suppressed, as expected in equilibrium, due to 
the instability of the system. Such fluctuations will drive 
the system to another point of stability of the phase dia- 
gram. The evolution in time of those fluctuations are at 
the heart of our discussion. 

For first-order phase transitions, there can be two kinds 
of instabilities that dominate the dynamics of the phase 
conversion [32| . If a homogeneous system is brought into 
instability close enough to the coexistence line of the 
phase diagram, its dynamics will be dominated by large- 
amplitude, small-ranged fluctuations. In these cases, 
large amplitudes are necessary for the development of 
the phase transition once the system is in a metastable 
equilibrium. Those are usually referred to as bubbles (or 
droplets) and the process that creates them is called nu- 
cleation. In the other case, if the external perturbation 
is big enough and the system finds itself far from the co- 
existence line, the dominant fluctuations will have small 
amplitudes and large wavelengths, the process of phase 
conversion is called spinodal decomposition. In this work, 
we focus on thermal nucleation of quark matter as nu- 
clear matter is compressed in a stellar collapse, leaving a 
discussion on a possible role for spinodal decomposition 
and quantum nucleation to the final section. 

A standard, field-theoretical approach for thermal nu- 
cleation in one-component metastable systems was devel- 
oped by Langcr in the late sixties [s^ . In this formalism, 
a key quantity for the calculation of the rate of nucleation 
is the coarse-grained free energy functional 



m = j 



VW)] 



(10) 



where 0(r) is the order parameter of the phase transition 
at a given point r of space. By assumption, the "po- 
tential" V{4i) has a global (true) minimum at 0t and 
a local (false) one at (/i/. At a given baryon chemi- 
cal potential fi of the metastable phase, the difference 
AF = V{(t)t) — V{4)f) is identified with the pressure dif- 
ference between the stable and the metastable phases, 
with opposite sign: Ay = —Ap{fj,) = pt — Pf, where pt 
{pf) is the pressure for the true (false) phase at baryon 
chemical potential fi. 

The field equation for (f)(r) is given by a minimum of 
the functional F. One can easily think of three static so- 
lutions. Two of them are the trivial ones given by homo- 
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gcneous field configurations with c/)(r) = (j)t or (j){r) = (/)/■ 
The third is a spherically symmetric bubblelike solution 
that has as boundary conditions 



= 0) = (j>t, 
(r go) = (j>f. 



(11) 



Roughly speaking, this means that the stable phase is 
found deep in the bubble and the mctastablc one is found 
away from it. Somewhere in-between, the order parame- 
ter must change from its central value (pt to cf>f at r oo. 
The relatively thin region which marks the border be- 
tween "inside" (0 = 4't) and "outside" (0 = (j)f) the 
bubble is called the bubble wall. 

Exactly at the coexistence line, one can prepare one 
(infinite) system with the two homogeneous phases in 
equal proportions divided by a plane wall with a small 
width. This configuration is static, once no phase is fa- 
vored. Further, each phase occupies a semi-infinite vol- 
ume. If the system is slightly pushed into metastability, 
the static solution for c/)(r) is a bubble with a very large 
radius and still a small wall width. This is the starting 
point for the thin-wall approximation: the free energy 
PU)) of the system of volume {4:n/3)L^ [L oo) is deter- 
mined by the outcome of a competition between a surface 
energy term, which is positive and comes from |V(/)p in 
pO)) . and a bulk term, which is negative and corresponds 
to the potential y, or to the pressure difference between 
the phases. Notice that, within this approximation, 0(r) 
is constant, except over the (thin) wall of the bubble, 
and so V{(j)) is also essentially constant both inside and 
outside the bubble. This means that the free energy for 
the bubble configuration of radius R in the thin-wall ap- 
proximation of Eq. (jlO[) is given by 

FhubUciR) - 4^i?V - — (i3 _ r3)p^ _ _i?3p^ (12) 

whereas the homogeneous metastable configuration has 

|V(/)|2 ^ and 
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In Eq. (|12p. we introduced the surface tension cr, which 
is merely the energy per unit area of the bubble wall. As 
will be clear below, it is a key physical quantity in our 
analysis. 

According to the standard theory [3^, the nucleation 
rate has as its main ingredient the free energy shift when 
a bubble is created from fluctuations in the homogeneous 
metastable phase. From to Eqs. ([T^ and we have 



AF{R) — ^bubblc(^) — ^'metastable 

= Attract R^{Ap), 



(14) 



where Ap = Pt—Pf > 0. Here, the pressures in each of the 
phases are calculated for the same value of Heff- Notice 
that this implies different baryon chemical potentials and 
densities for each phase, due to the conditions ([S])-®. 



Bubble configurations of given radii R arise from the 
homogeneous metastable phase due to thermal fluctu- 
ations, and each of those has an associated value of 
AF{R). From Eq. ((HI), we can see that AF{R) has a 
maximum at the critical radius Rc = 2a/ Ap. The equa- 
tions of motion show that any bubble with R < Rc will 
shrink and disappear whereas any bubble with R > Rc 
will grow, as a consequence of the competition between 
the positive surface energy and the negative bulk energy. 
Hence, the critical bubbles are the smallest bubbles that 
can start to drive the phase conversion dynamics. To 
give a quantitative meaning to the process of nucleation, 
one can calculate the rate F of critical bubbles created by 
fluctuationsper unit volume, per unit time. In Langer's 
formalism 1331: 
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where the prefactor 'Pg is usually factorized into two 
parts: a statistical prefactor, which measures the rate of 
successful creation of a critical bubble by thermal fluc- 
tuations, and a dynamical prefactor, which measures the 
early growth rate of the bubble. As customary, we adopt 
VqI^t: = T*, which corresponds to an overestimate of 
the actual prefactor. (For an exact calculation of Vq see, 
e.g., Ref. [SJ].) This constitutes one of our main reasons 
to interpret our results as providing an overestimate for 
the nucleation rate. It goes in line with the thin-wall ap- 
proximation, which is also known to overestimate F when 
compared to the exact (numerical) result [1^ . Although 
this overestimate can lead to an overall factor of ^ 10^ 
or even higher [sj] , the qualitative aspects of the results 
shown in the next section can be barely changed. And, 
since we are concerned with providing underestimates for 
thermal nucleation time scales under core collapse super- 
novae typical conditions, these details are not relevant. 
Our final formula for the nucleation rate reads 



F = exp 



167r 



3 (Ap)2T 



(16) 



where we used Eqs. and Rc ~ 2a/ Ap. Notice that 
the influence of the equation of state is present through 
Ap. Also, there is a remarkably strong dependence of F 
on the surface tension cr, which will be determinant for 
the nucleation time scale. 

It is convenient to introduce the nucleation time t, 
defined as the time it takes for the nucleation of one single 
critical bubble inside a volume of Ifcm'^ inside the core of 
the protoneutron star: 



1 



1 



Ikm^ J F 



(17) 



This is the time scale to be compared with the duration 
of the early post-bounce phase of a supernova event, few 
hundreds of milliseconds, during which it has been shown 
that quark matter formation could trigger the explosion 
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[6|. With this definition, we assume that temperature 
and density are constant within this central volume of 
Ikm^. Of course, this also goes in the direction of un- 
derestimating this time scale. In a more realistic calcu- 
lation, one should first compute the pressure and density 
profiles using the TOV equations, then calculate the local 
value of r, and finally integrate over the region contain- 
ing mctastable matter. However, the density profiles arc 
almost flat within the central kilometers of the star, thus 
making our assumption quite reasonable. 

Finally, we note that we calculate the time of produc- 
tion of one single critical bubble, which has a typical 
size of some fermi. We do not study, in this paper, the 
growth regime of the quark front. Once again, this leads 
to an underestimate: in comparing r with the bounce 
time scale as a criterion for the formation of a quark core, 
we tacitly assume that the quark matter bubble becomes 
macroscopic almost instantaneously, an obviously artifi- 
cial simplification. 



C. The role of statistical fluctuations 

As discussed in the Introduction, a possible alterna- 
tive approach to model the phase transition to quark 
matter assumes isospin (and strangeness) conservation 
during the phase transition (lH . [20I . [2ll |. This implies 
a transition from hadronic matter in chemical equilib- 
rium to an intermediate quark matter phase Q* in which 
quarks are not in chemical equilibrium. The chemical 
potentials of quarks are calculated, indeed, by assuming 
that the fractions of different quark flavors are the same, 
both in the hadronic and in the quark phase. This as- 
sumption is based on the argument that the time scale of 
the phase transition is regulated by QCD, thus typically 
of the order of 10"^'^ s, and much faster than the weak 
interaction time scale (weak interactions would produce 
quark matter in chemical equilibrium Q'"'^ only after the 
phase transition is completed). 

However, as noticed in Refs. (H,!!^, statistical fluctu- 
ations of the number densities of quarks can have impor- 
tant effects on nucleation. To estimate these effects one 
can calculate the phase transition point by considering 
the two equations of state in chemical equilibrium. At 
the critical density one compares the average numbers 
^u'd ^^"-^ down quarks within Q^'^ and within the 

nuclear phase N* in a flxed volume V. If these num- 
bers are different only by 3y^7V* ^, it means that most 

probably fluctuations would drive the transition directly 
to Q^"?. 

The volume V in which we consider fluctuations cor- 
responds to the volume of a drop of the new phase with 
critical radius Rcrit = 2a/ AP. Again, the surface ten- 
sion a is the crucial quantity that determines whether the 
phase transition occurs via the intermediate phase Q* or 
directly to Q'^'^. Taking into account the uncertainties on 
the value of a we estimate the critical radii to be of order 



^ 6 fm. Using the simple calculation explained above, we 
obtain that statistical fluctuations are efficient for radii 
of the order of 2 — 4 fm and thus of the same order of 
magnitude as the critical radii. So, we assume, as in Ref. 
[3], that the first drop is nucleated already in the 
phase. 

Now that we have all the ingredients for the calculation 
of the nucleation times, we can proceed to our results in 
the next section. 



III. RESULTS AND DISCUSSIONS 

In order to evaluate if the time scale r for the nucle- 
ation of quark matter is compatible with the bounce time 
scale tb, we underestimate the time r for the formation 
of a critical bubble as a function of density under var- 
ious conditions of temperature, as well as for different 
equations of state and values of surface tension, given 
the scenario described in the previous section. 

Recent supernova simulations indicate that the cen- 
tral density of a protoneutron star can be as high as 2no 
during the bounce, and this value will serve as a cutoff 
density in our analysis. Still in the spirit of underes- 
timating T (or, cquivalently, overestimating the nucle- 
ation rate F), we consider that nucleation is effective if 
T < Tb < 100 ms for some n < 2no. 



A. Nucleation times for nonstrange matter 

As our first case, we consider the transition from beta- 
stable nuclear matter to beta-stable quark matter com- 
posed of u and d quarks, plus electrons and electron neu- 
trinos, with a fixed lepton fraction Y^. Later on, we will 
discuss the case of a transition from nuclear matter to 
u-d-s quark matter (both lepton-rich and in beta equilib- 
rium). 

We start our analysis by the case of a low-density tran- 
sition: ric = 1.5no. We assume a lepton fraction = 0.4 
and consider two values for the temperature, representing 
a "minimum" and a "maximum" value that can be ex- 
pected during the bounce, and two values for the surface 
tension. 

In Fig. 1 one can see the behavior of the nucleation 
time of a single critical bubble (as defined in the previous 
section) versus the density, in units of uq. 

As expected, the nucleation time t has an extremely 
strong dependence on both density (notice the logarith- 
mic scale for r) and on the surface tension, a feature that 
can also be seen in Fig. 2. For low values of a nucleation 
becomes feasible at relatively low densities, although such 
densities increase steadily as the surface tension rises. 
However, if the (basically unknown) surface tension is 
larger, the density for r ^ 100 ms may be higher than 
our 2nQ cutoff, and nucleation should not be an efficient 
mechanism for phase conversion. In this sense, we can 
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FIG. 1: Nucleation time as a function of baryon density for u- 
d quark matter {ric = 1.5no). Tlie liorizontal line corresponds 
to r = 100ms. 
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FIG. 3: Lines of constant nucleation time (r = 100 ms) for 
Uc/no = 1.2, 1.5 , 2.0 with T = 20 MeV, c = and Yz, = 0.4. 
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FIG. 2: Contour lines of constant nucleation time (contour 
lines) as a function of density and surface tension for u-d 
quark matter. Here, iic ~ 1.5no, Yl = 0.4, and T = 20 MeV. 
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FIG. 4: Lines of constant nucleation time (r = 100ms) for 
noninteracting u-d quarks, for temperatures between 10 and 
25 MeV {ric = 1.5no, Yl = 0.4). 



expect that if the nuclear-quark matter transition occurs 
in this scenario the surface tension must be quite small. 

We can also investigate how the choice of the critical 
density can affect the nucleation time (Fig. 3). In any 
case, a surface tension larger than roughly 15 MeV/fm^ 
seems to be sufficient to prevent thermal nucleation (for 
T = 20 MeV or lower). 

In Fig. 4, we show the role played by the tempera- 
ture in the process of nucleation. We may notice that 
the precise value of the temperature does not affect the 
nucleation times as strongly as the surface tension does, 
as long as it is kept in the range expected during the 
early post-bounce phase at the protoneutron star core, 
i.e., roughly from 10 to 25 MeV. 

Although the exact numbers should not be taken at 
face value, given the uncertainties involved, we believe 



that the order of magnitude of the actual limiting value 
of the surface tension for r — 100 ms is correct. 



B. Nucleation of strange matter 

Although the core of a supernova progenitor star be- 
fore its collapse does not contain any strangeness, the 
energy density achieved during and right after bounce 
allows for the presence of a small amount of hyperons in 
the hadronic phase [37| . Such particles do not contribute 
significantly to the pressure or to the energy density, but 
density fluctuations of such hadrons may induce the for- 
mation of bubbles of strange quark matter. 

The introduction of strange quarks makes the EoS 
stiffer, i.e. for a given baryon chemical potential fi the 
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corresponding pressure becomes higher. Once the nu- 
clear EoS is the same, Ap will be higher for a given value 
of fi, and therefore the nucleation rate will also be higher. 
Figure 5 shows a comparison between the transition from 
nuclear matter to either u-d or u-d-s quark matter for two 
values of the lepton fraction Yl- We can notice that a 
decrease in Yl increases the efficiency of thermal nucle- 
ation. This is expected because dclcptonization not only 
renders nuclear matter less stable but it also heats up the 
stellar core (the former effect, however, is not accounted 
for in this work). 




FIG. 5: Lines of constant nucleation time (r = 100 ms) for the 
transition from nuclear matter to u-d or u-d-s quark matter, 
with Yl = 0.3, 0.4 (n^ = 1.5no, T = 20 MeV, c = and 
ms = 0). 



Up to now, we have only considered noninteracting 
quarks. Results from two-loop pcrturbative three-flavor 
QCD at finite density [25| show that strong interactions 
can be effectively accounted for in the equation of state 
by introducing a factor c < 1 [2^, as in Eq. [T] This fac- 
tor makes the quark EoS softer in the pressure-chemical 
potential plane [see Eq. dT])] and, therefore, Ap should be 
smaller, making the nucleation time r larger for a given 
density, according to Eqs. (fT6|) and (fT7|). As an explicit 
example, we compare the c = case with c = 0.3, in 
the case of Uc = 1.5rio, as displayed in Fig. 6, where we 
also show the influence of strange quark mass . Notice 
that the introduction of interactions via the parameter c 
drastically increases the nucleation time, so that only for 
a low value of the surface tension, e.g. a ^ 10 MeV/fm^ 
(for Yl = 0.4), nucleation can be efficient if the density 
reaches a value close to 2nQ. ^ 



c=0 


m 


=0 


--c=0 


m 


=100MeV 


■ - ■ c=0.3 


m 


=0 


■ ■ ■ ■ c=0.3 


m 


=100MeV 



.2 20- 

a 

<o 
H 

^ 10- 
00 
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FIG. 6: Lines of constant nucleation time (r = 100 ms) for 
c = and c — 0.3, and for nis = and nis — 100 MeV 
{uc = 1.5no, T = 20 MeV and Yl = 0.4). 



C. Stellar structure of selected equations of state 

In order to check if the equations of state we used are 
compatible with observed pulsar data, we calculate their 
associated mass-radius diagram using the TOV equa- 
tions. In Table 1, we report the quark model parameters, 
for strange matter, Yl = 0.4 and different choices of the 
critical density and the resulting maximum mass for cold 
and beta-stable stars. 

Notice that only by including the effect of QCD per- 
turbative interactions is it possible to obtain masses for 
hybrid stars compatible with the recent observation of 
PSR J1903-F0327, M = (1.671 ± 0.008)Mo 



Uc/nQ 


c 


^1/4 (MeV) 




1.2 





159.22 


1.60 




0.3 


144.65 


1.90 


1.5 





161.77 


1.55 




0.3 


145.89 


1.87 


2.0 





166.64 


1.48 




0.3 


147.56 


1.83 



Table 1. Maximum masses of cold deleptonized com- 
pact stars for some of the EoS used (corresponding to 
the cases Yl = 0.4 and nis = 100 MeV considered for 
nucleation). 

D. Quantum nucleation and spinodal instability 



^ Of course, the surface tension should also be affected by loop 
corrections, and that could eventually reduce or even balance out 
this effect on the final nucleation dynamics. As becomes clear 
from all this analysis, a reliable estimate of the surface tension 
for cold dense matter is called for. 



There are other mechanisms which can compete with 
thermal nucleation in driving the phase transition: quan- 
tum nucleation and spinodal decomposition. Unfortu- 
nately, our treatment is blind to the spinodal instability: 
an effective potential is indeed needed to take it into ac- 
count effectively (see Ref. [l^). In any case, this is a 
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process that will be relevant only if the system is taken 
into values of density that are high enough to flatten out 
the activation barrier, so that there is no more extra cost 
to create a region of the true vacuum inside the initially 
homogeneous false vacuum configuration by thermal fluc- 
tuations. If that is possible, the phase conversion process 
will be rather explosive, a possibility we will consider in 
a forthcoming publication. 

On the other hand, we can easily estimate the con- 
tribution of quantum nucleation by using the formalism 
considered in Ref. [l^l , i.e. a WKB treatment of the tun- 
neling through a barrier of an effective potential similar 
in form to Eq. ([T4)) . We find that only for temperatures 
smaller than ^ 5 MeV the quantum nucleation rate is 
comparable or larger than the thermal nucleation rate. 
Since we are considering supernova matter, with T > 10 
MeV, thermal nucleation is by far the dominant mecha- 
nism for the production of the first drop of quark matter. 

IV. CONCLUSIONS 

We investigated the possibility of the formation of 
quark matter in supernova matter, i.e. for temperatures 
of the order of a few tens of MeV and in the presence 
of trapped neutrinos, assuming that the corresponding 
critical density does not exceed 2no- We argued that 
thermal nucleation of quark phase droplets is eventually 
the dominant mechanism for the formation of the new 
phase and that, due to fiuctuations in the number densi- 
ties of quarks, the phase transition involves directly the 
beta equilibrated quark phase. We have calculated the 
nucleation rate for different values of the free parameters. 
The surface tension, as expected, is the physical quantity 
which mainly controls the nucleation process. 

Among the different equations of state and conditions 
at bounce we have tested, the choice T = 20 MeV and 
Yl = 0.4 is the most likely to occur. Within this choice of 
physical conditions, if the phase transition involves only 



noninteracting up and down quarks, a value of a smaller 
than ~ 15 MeV/fm^ should be required for nucleation 
to be efficient. Such a low value of a is compatible with 
lattice QCD calculations, although they are obtained at 
large temperatures and small densities [38| . On the other 
hand, phenomenological estimates at large density and 
zero temperature indicate larger values of a (soj . There- 
fore, we consider this scenario to be unlikely. 

If strange quarks are produced during the phase transi- 
tion, we conclude that, if a is smaller than ~ 10 MeV/fm^ 
(for nis — lOOMeV and c = 0.3), the nucleation time for 
the first drop of quark matter is sufficiently small and 
the appearance of quark matter can indeed strongly af- 
fect the supernova evolution as shown in Q. The cold 
hybrid stars obtained after deleptonization and cooling 
have, if perturbative QCD corrections are included in 
the equation of state, maximum masses compatible with 
recent pulsar observations. 

Note added- After finishing this work a paper has been 
published which discusses similar issues [401. In. that 
work, the authors also conclude that nucleation is pos- 
sible in protoneutron star matter. The main difference 
between our work and Ref. [i^l is that in the latter the 
thermal nucleation of quark matter is studied within hot 
and deleptonized hadronic matter. 
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